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Abstract 

It is known that the sphere packing problem and the number variance problem (closely related 
to an optimization problem in number theory) can be posed as energy minimizations associated 
with an infinite number of point particles in d-dimensional Euclidean space M d interacting via 
certain repulsive pair potentials. We reformulate the covering and quantizer problems as the 
determination of the ground states of interacting particles in M. d that generally involve single- 
body, two-body, three-body, and higher-body interactions. This is done by linking the covering 
and quantizer problems to certain optimization problems involving the "void" nearest-neighbor 
functions that arise in the theory of random media and statistical mechanics. These reformulations, 
which again exemplifies the deep interplay between geometry and physics, allow one now to employ 
theoretical and numerical optimization techniques to analyze and solve these energy minimization 
problems. The covering and quantizer problems have relevance in numerous applications, including 
wireless communication network layouts, the search of high-dimensional data parameter spaces, 
stereotactic radiation therapy, data compression, digital communications, meshing of space for 
numerical analysis, and coding and cryptography, among other examples. In the first three space 
dimensions, the best known solutions of the sphere packing and number variance problems (or 
their "dual" solutions) are directly related to those of the covering and quantizer problems, but 
such relationships may or may not exist for d > 4, depending on the peculiarities of the dimensions 
involved. Our reformulation sheds light on the reasons for these similarities and differences. We also 
show that disordered saturated sphere packings provide relatively thin (economical) coverings and 
may yield thinner coverings than the best known lattice coverings in sufficiently large dimensions. 
In the case of the quantizer problem, we derive improved upper bounds on the quantizer error 
using sphere-packing solutions, which are generally substantially sharper than an existing upper 
bound in low to moderately large dimensions. We also demonstrate that disordered saturated 
sphere packings yield relatively good quantizers. Finally, we remark on possible applications of our 
results for the detection of gravitational waves. 

PACS numbers: 05.20.-y, 89.20.Ff, 89.70.-a, 89.75.Kd 
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I. INTRODUCTION 



There are certain scientific problems that provide deep connections between many differ- 
ent scientific fields. The study of the low-energy states of classical interacting many-particle 
systems is an exemplar of a class of such problems because of its manifest importance in 
physics, materials science, communication theory, cryptography, mathematics and computer 
science. Such many-particle systems have been used with great success to model liquids, 
glasses and crystals when quantum effects are negligible l],!^. The total potential energy 
$Ar(r Ar ) of N identical particles with positions = r 1; r 2 , . . . , rjy in some large volume in 
d- dimensional Euclidean space M. d can be resolved into separate one-body, two-body, . . . , 
iV-body contributions: 

N N N 

® N (r N ) =^w 1 (r i ) + ^« 2 (ri,r J -)+ u 3 {r i ,r j ,r k ) + --- + u N (r N ), (1) 

i=l i<j i<j<k 

where u n represents the intrinsic n-body interaction in excess to the interaction energy for 
n — 1 particles. To make the statistical-mechanical problem more tractable, the exact many- 
body potential is usually replaced by a mathematically simpler form. For example, in 
the absence of an external field (i.e., u\ =0), often one assumes pairwise additivity, i.e., 



N 



(2) 



Pairwise additivity is exact for hard-sphere systems and frequently has served to approximate 
accurately the interactions in simple liquids, such as the well-known Lennard- Jones pair 
potential [l|, and in more complex systems where the pair potential u<i in ([2]) can be regarded 
to be an effective pair interaction 3J. 

An outstanding problem in classical statistical mechanics is the determination of the 
ground states of <&n(t n ), which are those configurations that globally minimize $>n(t n ) and 
hence are the states that exist at absolute zero temperature. While classical ground states 
are readily produced by slowly freezing liquids in experiments and computer simulations, 
our theoretical understanding of them is far from complete 0, 0]. Virtually all theoreti- 
cal/computational ground-state studies of many-particle systems have been conducted for 
pairwise additive potentials [l|, Q, E ll]. Often the ground states of short-range pairwise 



interactions are crystal structures in low dimensions [l|, |4|, [6 



10 



121 ]. but long-range inter- 



actions exist that can suppress any kind of symmetry leading to disordered ground states 



in low dimensions [kJQ]. Moreover, in sufficiently high dimensions, it has been suggested 
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that even short-ranged pairwise interactions possess disordered ground states 

Ground states of purely repulsive pair interactions have profound connections not only to 
low-temperature states of matter but to problems in pure mathematics, including discrete 
geometry and number theory 1Q,Q)Q|) information theory, and computer science. As will 
be explained further below, it is known that the sphere packing problem and the number 
variance problem (closely related to an optimization problem in number theory) can be 
posed as energy minimizations associated with an infinite number of point particles in d- 
dimensional Euclidean space M. d interacting via certain repulsive pair potentials. Both of 
these problems can be interpreted to be optimization problems involving point processes, 
which can then be recast as energy minimizations. A point process in M. d is a distribution 
of an an infinite number of points in M. d at number density p (number of points per unit 
volume) with configuration r 1; r 2 , . . .; see Ref. 16[ for a precise mathematical definition. 

A packing of congruent nonoverlapping spheres is a special point process in which there 
is a minimal pair separation distance, equal to the sphere diameter. The sphere packing 
problem seeks to determine the densest arrangement (s) of congruent, nonoverlapping d- 



dimensional spheres in Euclidean space 



18 



201 ] . Although it is simple to state, it is a 



notoriously difficult problem to solve rigorously. Indeed, Kepler's four-century-old conjec- 
ture, which states that the face-centered-cubic lattice in M. 3 is maximally dense, was only 
recently proved 2l|. For d > 4, the packing problem remains unsolved 18|, |20|, |22| . It is well 
known that the sphere packing problem can be posed as an energy minimization problem 
involving pairwise interactions between points in M. d (e.g., inverse power-law functions in 
which the exponent tends to infinity); see Ref. 10|] and references therein. 

Problems concerning the properties and quantification of density fluctuations in many- 
particle systems continue to provide many theoretical chal leng es. Of particular interest are 
density fluctuations that occur on some local length scale 23[. It has been shown that the 
minimal number variance associated with points (e.g., centroids of atomic or molecular sys- 
tems) contained within some "window" can also be formulated as a ground-state problem 
involving bounded repulsive pair interactions with compact support 23|, |24|. For spherical 
windows in the large-radius limit, the best known solutions in M d are usually point configu- 
rations that are "duals" (in the sense discussed later in the paper) to the best known sphere 
packings in M. d . 
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The focus of this paper is on two other optimization problems involving point processes 
in M. d : the covering and quantizer problems. Roughly speaking, the covering problem asks 
for the point configuration that minimizes the radius of overlapping spheres circumscribed 
around each of the points required to cover R d . The covering problem has applications in 
wireless communication network layouts 25j, the search of high-dimensional data parameter 
spaces (e.g., search templates for gravitational waves) 26], and stereotactic radiation therapy 



27| . The quantizer problem is concerned with finding the point configuration in M. d that 



minimizes a "distance error" associated with a randomly placed point and the nearest point 
of the point process. It has applications in computer science (e.g., data compression) 18], 



digital communications 



18| . coding and cryptography [28| . and optimal meshing of space 



or numerical applications (e.g., quadrature and discretizing partial differential equations) 



291 ] . Heretofore, the covering and quantizers problems were not known to correspond to any 



ground-state problems. 

We reformulate the covering and quantizer problems as the determination of the ground 
states of interacting particles in M. d that generally involve single-body, two-body, three-body, 
and higher-body interactions. This is done by linking the covering and quantizer problems 
to certain optimization problems involving the "void" nearest-neighbor functions that arise 
in the theory of random media and statistical mechanics {2], 3o|, 31 1. These reformulations, 
which again exemplifies the deep interplay between geometry and physics, enable one to em- 
ploy theoretical and numerical optimization techniques to solve these energy minimization 
problems. We find that disordered saturated sphere packings (roughly, packings in which 
no space exists to add an additional sphere) provide relatively thin (i.e., economical) cov- 
erings and may yield thinner coverings than the best known lattice coverings in sufficiently 
large dimensions. In the case of the quantizer problem, we derive improved upper bounds 
on the quantizer error that utilize sphere-packing solutions. These improved bounds are 
generally substantially sharper than an existing upper bound in low to moderately large 
dimensions. We also demonstrate that disordered saturated sphere packings yield relatively 
good quantizers. Our reformulation helps to explain why the known solutions of quantizer 
and covering problems are identical in the first three space dimensions and why they can 
be different for d > 4. In the first three space dimensions, the best known solutions of the 
sphere packing and number variance problems are directly related to those of the covering 
and quantizer problems, but such relationships may or may not exist for d > 4, depending 
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on the peculiarities of the dimensions involved. 

We begin by summarizing basic definitions and concepts in Sec. [XTJ Because of the 
connections between the sphere-packing, number- variance, covering and quantizer problems, 
in Sec. IHIl we formally define each of these problems, summarize key developments, and 
compare the best known solutions for each of them in selected dimensions. This includes 
calculations obtained for the best known number-variance solutions for d = 12, 16 and 24. 
We then define in Sec. [IV] the void nearest-neighbor functions and represent them in terms of 
series involving certain integrals over the n-particle correlation functions, which statistically 
characterize an ensemble of interacting points. The special case of a single realization of 
the point distribution follows from this ensemble formulation, which reveals that quantizer 
and covering problems can be expressed as ground-state solutions of many-body interactions 
of the general form ([1]). Section |V] specifically gives these ground-state reformulations and 
shows how some known solutions in low dimensions can be explicitly recovered using the void 
nearest-neighbor functions. In Sec. IVIl we show that disordered saturated sphere packings 
provide relatively thin coverings and may yield thinner coverings than the best known lattice 
coverings in sufficiently large dimensions. In Sec. IVIIl we derive improved upper bounds on 
the quantizer error that utilize sphere-packing solutions. We also show that disordered 
saturated sphere packings yield relatively good quantizers. Finally, in Sec. IVHH we make 
concluding remarks and comment on the application of the quantizer problem to the search 
for gravitational waves. 

II. DEFINITIONS AND PRELIMINARIES 

For a statistically homogeneous point process in M. d at number density p (number of 
points per unit volume), the quantity p n g n (ri, r 2 , . . . , r n ) is proportional to the probability 
density for simultaneously finding n sphere centers at locations ri, ■ ■ ■ , r n in M. d [l|. With 
this convention, each n-particle correlation function g n approaches unity when all of the 
points become widely separated from one another. Statistical homogeneity implies that g n 
is translationally invariant and therefore only depends on the relative displacements of the 
positions with respect to some arbitrarily chosen origin of the system, i.e., 

g n = 9n{ r l2, 1"l3, • • • , Tin), (3) 
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where = rj — Tj. As we will see, statistically homogeneous point processes include as 
special cases periodic point distributions. 

The pair correlation function g 2 (r) is a particularly important quantity. If the point 
process is also rotationally invariant (statistically isotropic), then g 2 depends on the radial 
distance r = |r| only, i.e., <?2( r ) = 92(j')- Thus, it follows that the expected number of points 
Z(R) found in a sphere of radius R from a randomly chosen point of the point process, 
called the cumulative coordination function, is given by 



is the surface area of a rf-dimensional sphere of radius r. 

A lattice A in M. d is a subgroup consisting of the integer linear combinations of vectors that 
constitute a basis for JBL d and thus represents a special subset of point processes. In a lattice 
A, the space M. d can be geometrically divided into identical regions F called fundamental 
cells, each of which contains the just one point specified by the lattice vector 



where a^ are the basis vectors for the fundamental cell and rti spans all the integers for 
i — 1, 2, • • • d. We denote by vp the volume of the fundamental cell. In the physical sciences, 
a lattice is equivalent to a Bravais lattice. Unless otherwise stated, we will use the term 
lattice. Every lattice has a dual (or reciprocal) lattice A* in which the sites of the lattice 
are specified by the dual (reciprocal) lattice vector q • p = 2-Km, where m — ±1, ±2, ±3 • • • . 
The dual fundamental cell F* has volume vp* = (2n) d /vp. This implies that the number 
density p of A is related to the number density of the dual lattice A* via the expression 
pp* = l/(27r) d . A periodic point process is a more general notion than a lattice because 
it is is obtained by placing a fixed configuration of N points (where N > 1) within one 
fundamental cell of a lattice A, which is then periodically replicated. Thus, the point process 
is still periodic under translations by A, but the iV points can occur anywhere in the chosen 
fundamental cell. 




(4) 



where 




(5) 



p = mai + n 2 a 2 H h n d „ia d _i + n d & d , 



(6) 
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Common d- dimensional lattices include the hypercubic Z d , checkerboard Dd and root A d 
lattices, denned, respectively, by 

Z d = {(an, . . . , x d ) : Xi e Z} for d > 1 (7) 

Dd = {{xi-, ■ ■ ■ i Xd) € Z d : x\ + ■ ■ • + Xd even} for d > 3 (8) 

A d = {{x ,xi,. . .,x d ) E Z d+1 : x +x 1 + • •• + x d = 0} for d > 1 (9) 

where Z is the set of integers (. . . — 3, —2, —1, 0, 1, 2, 3 . . .) and x±, . . . , x d denote the compo- 
nents of a lattice vector of either Z d or D d and Xq, Xi, . . . ,x d denote a lattice vector of Ay. 
The d- dimensional lattices Z d , D* d and A* d are the corresponding dual lattices; see Ref. [18] 



for definitions. The dual lattice Z d is also a hypercubic lattice (even if the lattice spacing 
is 2tt times the lattice spacing of Z d ) and hence we say that the hypercubic lattice Z d is 



equivalent (similar) to its dual lattice Z^, i.e., 7L d = Z^. Following Conway and Sloane 181 ] . 
we say that two lattice are equivalent or similar if one becomes identical to the other by 
possibly a rotation, reflection and change of scale, for which we use the symbol =. In fact, 
the hypercubic lattice Z d is characterized by the stronger property of self-duality. A self- 
dual lattice is one with an identical dual lattice at density p = p* = l/(27r) d//2 , i.e., without 
any rotation, reflection, or change of scale [3^]. The A d and Dd lattices can be regarded to 
be d- dimensional generalizations of the face-centered-cubic (fee) lattice because this three- 
dimensional lattice is defined by A3 = D%. In one dimension, A\ = A* (equality meaning 
self-duality) are identical to the integer lattice Z 1 = Z. In two dimensions, A 2 = A* 2 defines 
the triangular lattice. In three dimensions, A\ = D% defines the body-centered-cubic (bec) 



IS 



331 ] is of special interest. In dimensions 



lattice. The d-dimensional laminated lattice 
8 and 24, A 8 = E 8 , where E s is the self-dual root lattice, and A 2 4 is the self-dual Leech 
lattice, are remarkably symmetric and believed to be the densest sphere packings in those 
dimensions [22]. Thus, E$ = Eg and A24 = A24. The laminated lattice Ai6, called the 
Barnes- Wall lattice, and the Coxeter-Todd lattice K 12 are thought to be the densest lattice 
packings in sixteen and twelve dimensions, respectively. 

Note that for a single periodic point configuration at number density p, the radial pair 
correlation function can be written as 

00 „ 

*< r >=£^ 4(r " r,) ' (10) 



s 



where Z\ is the coordination number at radial distance (number of points that are exactly 
at a distance r = from a point of the point process) such that r^+i > r% and 5(r) is a radial 
Dirac delta function. For cases in which each point is equivalent to any other, which includes 
all lattices and some periodic point processes, the coordination numbers Zj are integers. For 
point processes for which the point are generally inequivalent, Zj should be interpreted as 
the expected coordination number and hence will generally be a non-integer. Substitution of 
f fTOj) into @ gives the coordination function for such a periodic configuration as 

M 

Z(R) = J2 Z " ( U ) 

t=i 

where M is the smallest integer for which tm+i > R- 

Consider any discrete set of points with position vectors Y = {ri, r 2 , . . .} in M. d . Associ- 
ated with each point rj G Y is its Voronoi cell, V(rj), which is defined to be the region of 
space nearer to the point at than to any other point Yj in the set, i.e., 

Vfo) = {x : |x - r<| < |x - t s \ for all r,- G Y}. (12) 

The Voronoi cells are convex polyhedra whose interiors are disjoint, but share common faces, 
and therefore the union of all of the polyhedra is the whole of M. d . This partition of space is 
called the Voronoi tessellation. The vertices of the Voronoi polyhedra are the points whose 
distance from the points Y is a local maximum. While the Voronoi polyhedra of a lattice are 
congruent to one another, the Voronoi polyhedra of a non-Bravais lattice are not identical to 
one another. A hole in a lattice is a point in M. d whose distance to the nearest lattice point is 
a local maximum. A deep hole is one whose distance to a lattice point is a global maximum. 
The distance 1Z C to the deepest hole of a lattice is the covering radius and is equal to the 
circumradius of the associated Voronoi cell (the radius of the smallest circumscribed sphere). 

In the case of the d- dimensional simple cubic (hypercubic) lattice Z d , the Voronoi cell is 
a hypercube and there is only one type of hole with covering radius 1Z C = y/d/2, assuming 
unit number density p— 1. Figure [1] shows the Voronoi cell in the case d = 3 as well as the 
corresponding Voronoi cells for the three-dimensional body-centered cubic and face-centered 
cubic lattices. Note that the truncated octahedron is the most spherically symmetric of the 



34|. 



three Voronoi cells shown in Fig. [T] 

A sphere packing P in d- dimensional Euclidean space M. d is a collection of d- dimensional 
nonoverlapping congruent spheres. The packing density or, simply, density 4>{P) of a sphere 



FIG. 1: (Color online) Voronoi cells in M 3 for simple cubic, body-centered cubic, and face-centered 
cubic lattices are the cube (left), truncated octahedron (middle), and rhombic dodecahedron (right), 



as adapted from Ref. 



2J. The truncated octahedron is composed of six square and eight regular 



hexagonal faces. The rhombic dodecahedron is composed of twelve rhombus-shaped faces. 

packing is the fraction of space M. d covered by the spheres. For spheres of diameter D and 
number density p, the density is given by 

<P = p Vl (D/2), (13) 

where 

<" w = WTW) Rd < 14) 

is the volume of a d- dimensional sphere of radius R. 

A packing is saturated if there is no space available to add another sphere without over- 
lapping the existing particles. We denote the packing density of such a packing by <p s . A 
lattice packing Pl is one in which the centers of nonoverlapping spheres are located at the 
points of A. Thus, the density of a lattice packing <ft L consisting of spheres of diameter D is 
given by 

? = (15) 
where vp is the volume of a fundamental cell. A periodic packing of congruent spheres is 
obtained by placing a fixed configuration of N sphere centers (where N > 1) within one 
fundamental cell of a lattice A, which is then periodically replicated without overlaps. The 
packing density of a periodic packing is given by 

, Nv 1 (D/2) 



where p = N/vp is the number density. 



pvx(D/2), (16) 
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III. PROBLEM STATEMENTS AND BACKGROUND 



A. Sphere Packing Problem 

The sphere packing problem seeks to answer the following question: Among all packings 
of congruent spheres, what is the maximal packing density max , i.e., largest fraction of R d 
covered by the spheres, and what are the corresponding arrangements of the spheres 
More precisely, the maximal density is defined by 



sup <f>(P), 



(17) 

where the supremum is taken over all packings in M d . The sphere packing problem is of great 
fundamental and practical interest, and arises in a variety of contexts, including classical 
ground states of matter in low dimensions 6J, the famous Kepler conjecture for d = 3 



error-correcting codes 



18 



361 ] and spherical codes 



211, 



2l|. For 



The optimal solutions are known only for the first three space dimensions 
4 < d < 9, the densest known packings are Bravais lattice packings [rsl ]. For example, 
the "checkerboard" lattice D^, which is the d- dimensional generalization of the fee lattice 
(densest packing in K 3 ), is believed to be optimal in M 4 and M 5 . The remarkably symmetric 
self-dual E$ and Leech lattices in M 8 and M 24 , respectively, are most likely the densest 
packings in these dimensions [22|. Table Q] summarizes the densest known packings in selected 
dimensions. 



For large (L 



on 



16 



the best that one can do theoretically is to devise upper and lower bounds 



18 



[20 ] . For example, Minkowski [37j proved that the maximal density 
among all Bravais lattice packings for d > 2 satisfies the lower bound 

((d) 



'max — 2^-1 



max 



(18) 



where ((d) = YlT=i k~ d is the Riemann zeta function. It is seen that for large values of d, 
the asymptotic behavior of the nonconstructive Minkowski lower bound is controlled by 2~ d . 
Note that the density of a saturated packing of congruent spheres in M. d for all d satisfies 

1 



> 



2 d ' 



(19) 



which has the same dominant exponential term as (1181) . This is a rather weak lower bound on 
the density of saturated packings because there exists a disordered but unsaturated packing 
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TABLE I: Best known solutions to the sphere packing problem in selected dimensions; see Conway 
and Sloane [18| for details. 



Dimension, d 


Packing 


Packing density, 


1 


A\ = Z 


1 


2 


a; = a 2 


7r/Vl2 = 0.906899 .. . 


3 


A 3 = D 3 


n/VT8 = 0.740480... 


4 


D A = Dt 


vr 2 /16 = 0.616850. . . 


5 


D 5 


2tt 2 /(30 v / 2) = 0.465257... 


6 


Eq 


3tt 3 /(144-v/3) = 0.372947... 


7 


E 7 


tt 3 /105 = 0.295297... 


8 


Es = Eg 


vr 4 /384 = 0.253669 . . . 


9 


A 9 


v/2tt 4 /945 = 0.145774... 


10 


P\0c 


tt 5 /3072 = 0.099615... 


12 


A max 
iv 12 


tt 6 /23040 = 0.041726 . . . 


16 


Aw 


tt 8 /645120 = 0.014708 . . . 


24 


A 24 = A24 


vr 24 /479001600 = 0.001929 . . . 



construction in M. d , known as the "ghost" RSA packing 15|, that achieves the density 2~ d in 
any dimension. We will employ these results in Sec. I VII Bl It is also known that there are 
saturated packings in M d with densities that exceed the scaling 2~ d [38[, as we will discuss 



in Sec. IIII CI In the large- dimensional limit, Kabatiansky and Levenshtein 
the maximal density is bounded from above according to the asymptotic upper bound 

CDm civ 



showed that 



(20) 



20.5990 d' 

We will employ sphere-packing solutions to obtain heretofore unattained results for both 
the covering and quantizer problems. In particular, we obtain coverings and quantizers 
utilizing disordered saturated packings in Sees. ED and I VII CI respectively. In Sec. IVHl we 
use the densest lattice packings to derive improved upper bounds on the quantizer error. 



B. Number Variance Problem, Hyperuniformity and Epstein Zeta Function 

We denote by o~ 2 (A) the variance in the number of points N(A) contained within a 
window A C IR d . The number variance c 2 (A) for a specific choice of A is necessarily a 
positive number and is generally related to the total correlation function h(r) = g 2 (r) — 1 



12 



for a translationally invariant point process 23j], where #2 is the pair correlation function 
defined in Sec. [TTJ In the special case of a spherical window of radius R in R d , the number 
variance is explicitly given by 



a 2 (R) =pv l {R) 



h(r)a(r; R) dr 



(21) 



a(r;R) is the dimensionless volume common to two spherical windows of radius R (in units 
of the volume of a spherical window of radius R, vi(R)) whose centers are separated by a 
distance r. We will call a(r; R) the scaled intersection volume, which will play an important 
role in this paper. The scaled intersection volume has the support [0, 2R], the range [0, 1], 



and the following alternative integral representation 



16|: 



/■cos" :L [r/(2_R)] 

a (r; R) = c(d) / sin d (#) d9, 

Jo 

where c(d) is the d- dimensional constant given by 

2r(l + d/2) 



c(d) 



(22) 



(23) 



7^/21^+ l)/2]- 

Torquato and Stillinger [16| found the following series representation of the scaled inter- 
section volume a(r; R) for r < 2R and for any d: 



a(r; it!) = 1 — c(d)x 



+ c(d) f :( _ ir (f-l)(^-3)---(^-2n + 3) 
1 ' (2n-l)[2.4.6...(2n-2)] 



(24) 



where x = r/(2R). For even dimensions, relation (|2~4"|) is an infinite series because it involves 
transcendental functions, but for odd dimensions, the series truncates such that a(r; R) is a 
univariate polynomial of degree d. For example, in two and three dimensions, respectively, 
the scaled intersection volumes are given by 

„2 \ 1/2" 



a(r; R) 



2 

Ti 



COS 



_j / r \ r 
\2~r) ~ 2R 



AR 2 



B(2R-r) (d = 2), 



a(r; R) 



where 



3 r 1 / r \ 
1 ~ 4R + 16 \R 



G(x) 



Q{2R-r) (d = 3), 




(25) 



(26) 



(27) 



is the Heaviside step function. Figure [2] provides plots of a(r; R) as a function of r for the first 
five space dimensions. For any dimension, a(r;R) is a monotonically decreasing function of 
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r/(2R) 



FIG. 2: The scaled intersection volume a(r; R) for spherical windows of radius R as a function of 
r for the first five space dimensions. The uppermost curve is for d = 1 and lowermost curve is for 
d = 5. 



r. At a fixed value of r in the interval (0, 2R), a(r; R) is a monotonically decreasing function 
of the dimension d. 

For large R, it has been proved that u 2 {R) cannot grow more slowly than 7i? d_1 , where 
7 is a positive constant |40|. We note that point processes (translationally invariant or not) 
for which <J 2 {R) grows more slowly than the window volume (i.e., as R d ) for large R are 
examples of hyperuniform (or superhomogeneous) point patterns {23 1. For hyperuniform 
point processes in which the number variance grows like the surface area of the window, one 
has 



a 2 {R) = 2 d r] 



B 



d-l 



O 



d-2 



R — y 00, 



where 



B 



77 c(d) 



2Vv 1 (V/2) 

is a dimensionless constant with c(d) defined by ( 123 



h(r)rdr, 



T] = pvtiV/2) 



(28) 



(29) 



(30) 



is a dimensionless density, and D represents some "microscopic" length scale, such as the 
minimum pair separation distance in a packing or the mean nearest-neighbor distance. This 
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class of hyperuniform point processes includes all periodic point patterns, quasicrystals that 
possess Bragg peaks, and disordered hyperuniform point patterns in which the pair correla- 
tion functions decay exponentially fast to unity 23]. 

It has been shown that finding the point process that minimizes the number variance 
<J 2 (R) is equivalent to finding the ground state of a certain repulsive pair potential with 



compact support 



231 ] . Specifically, by invoking a volume- average interpretation of the num- 



ber variance problem valid for a single realization of a point process, Torquato and Stillinger 
found [23J: 

d-i 



^ R) 2"/,Z?.v(/?)(^ 



where 



Bn(R) = § 



(31) 



(32) 



The asymptotic coefficient B defined by ff29l) for a hyperuniform point pattern is then related 
to B N (R) by the expression 



B 



lim — 

L->+oo L 



B N {R)dR. 



(33) 



These results imply that the asymptotic coefficient B obtained in ( 1291) involves an average 
over small-scale fluctuations in the number variance with length scale on the order of the 
mean separation between points 23|. In the special case of a (Bravais) lattice A, one can 
express the rescaled surface-area coefficient as follows: 



V 



i/d B 



7r (d-i)/22<*-i[r(i + d /2)] 1 - 1 /<2 



l+l/d 



E 



1 



\d+l : 



(34) 

V F q^O 11411 

where we recall that vp is the volume of the fundamental cell of the lattice A and q represents 
a lattice vector in the dual (or reciprocal) lattice A*. The rescaled coefficient r^l d B renders 



23) 



the result independent of the length scale in the lattice 

Finding the lattice that minimizes B is directly related to an outstanding problem in 
number theory, namely, finding the minima of the Epstein zeta function Z\(s) [l9[ defined 
by 

^a(s) - £ Tr^j, (35) 



where p is a lattice vector of the lattice A. Note that the dual of the lattice that minimizes the 
Epstein zeta function at s — (d+l)/2 among all lattices will minimize the scaled asymptotic 
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TABLE II: Best known solutions to the asymptotic number variance problem in selected dimen- 
sions. Values reported for d = 1 , 2 and 3 and d = 4-8 are taken from Torquato and Stillinger 23] 



and Zachary and Torquato [4l|], respectively. Values reported for d = 12,16 and 24 have been 
determined in the present work. 



Dimension, d 



1 

2 

3 

4 

5 

6 

7 

8 

12 

16 

24 



Structure 



A* = Z 
A* = A 2 
Al = D* 
D\ = D A 
Af 
E* 6 

Ea = Ea — Do 



K* 

-"■12 


= K n 


A* 


= Aie 


A* 


= A 24 



Scaled Coefficient, r) x l d B 



0.083333 
0.12709 
0.15560 
0.17488 
0.19069 
0.20221 
0.21037 
0.21746 
0.24344 
0.25629 
0.26775 



number- variance coefficient (|3"4"|) among lattices 23|,|4l|. Certain duality relations have been 
derived that establish rigorous upper bounds on the energies of such ground states and help 
to identify energy-minimizing lattices J5]. Because Za(s) is globally minimized for d = 1 



by the integer lattice 
lattice 



23| and is minimized for d = 2 among all lattices by the triangular 



421 ]. it has been conjectured that the Epstein zeta 



unction for s > is minimized 



among lattices by the — ly dease lattice packing U^. Sarnak and Strdmbergsson 
[191 ] have proved that the conjecture cannot be generally true, but for d = 4, 8 and 24, 
the densest lattice packing is a strict local minimum. Since as s — > +00, the minimizer of 
Epstein zeta function is the densest sphere packing in IR d for any d, it is likely that in the 
high-dimensional limit the minimizers of this function are non-lattices, namely, disordered 
sphere packings [161 ] . 

In Table HIl we tabulate the best known solutions to the asymptotic number variance 



problem in selected dimensions. Values for the first t 



iree space dimensions were given in Ref. 



231 ] and those for d = 4 — 8 were provided in Ref. 41] . The values reported for d = 12, 16 
and 24 were ascertained using efficient algorithms based on alternative number-theoretic 
representations of the Epstein zeta function Z\(s) 19j for the corresponding densest known 
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FIG. 3: Coverings of the plane with overlapping circles centered on the triangular lattice (left 
panel) and the square lattice (right panel). The triangular lattice A2 = A\ provides the best 
covering among all point processes at unit number density p = 1 with 9 = 2tt/(3\/3) = 1.2092 
This is to be compared to the square lattice 1? with 9 = vr/2 = 1.5708 .... 

lattice packings for s = (d + l)/2 and then using the duality relations connecting it to the 
asymptotic surface-area coefficient for the number variance. Appendix |A] provides details 
for these computations. 

C. Covering Problem 

Surround each of the points of a point process V in M, d by congruent overlapping spheres 
of radius R such that the spheres cover the space. The covering density 9 is denned as 
follows: 

9 = pv!{R), (36) 

where vi(R) is given by ( JT4l) . The covering problem asks for the arrangement of points with 
the least density 9. We define the covering radius TZ C for any configuration of points in M d to 
be the minimal radius of the overlapping spheres to cover M. d . Figure [3] shows two examples 
of coverings in the plane. 

The covering density associated with A* d at unit number density p = 1 is known exactly 



17 



for any dimension d 



9 = Vl (l)Vd+l 
For the hypercubic lattice 7L d at p = 1, 

= <ui(l) 



d(d + 2) 
12 (d + 1) 

d d ' 2 



Thus the ratio of the covering density for A* d to that of Z d is given by 



e(A* d ) Vd + i 



9{Z a 



3 d/2 



For large d, this ratio becomes 



0(A* d ) 



d + 2 
d+1 

de 



d/2 



(37) 



(38) 



(39) 



(40) 



9(Z d ) 3 d / 2 ' 

and thus we see that A* d provides exponentially thinner coverings than that of Z d in the 
large-d limit. We note that in this asymptotic limit, 



0{A* d ) 



e\V2 /e7r\i/2 

-J (_ J = 0.8652559792... (1.193016780...)- 



(41) 



Until recently, was the best known lattice covering in all dimensions d < 23. However, 
for most dimensions in the range 6 < 9 < 15, Schiirmann and Vallentin 45] have discovered 
other lattice coverings that are slightly thinner than those for A* d . Table II III provides the 
best known solutions to the covering problem in selected dimensions. 

Until the present work, there were no known explicit non-lattice constructions possessing 
covering densities smaller than those of the best lattice coverings in any dimension d Jja]. In 
Sec. IVII1 we provide evidence that certain disordered point patterns give thinner coverings 
than the best known lattice coverings beginning at d = 17. There is a fundamental difference 
between coverings associated with point patterns that have identical Voronoi cells (i.e., 
lattices and periodic point patterns in which each point is equivalent) and those point 
processes whose Voronoi cells are generally different (e.g., irregular point processes). This 
salient point is illustrated in Fig. |4] and explained in the corresponding caption. 

Rogers showed that (possibly nonlattice) coverings exist with 



9 < 5d + d\n(d) + dln(ln(d)). 



(42) 



for d > 2. This is a nonconstructive upper bound. This upper bound provides a substantially 
thinner covering density than that of the A* d lattice in the large-<i limit [cf. (j4"Tl) ]. but it is 
not known whether this bound becomes sharp in the large-d limit. 
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TABLE III: Best known solutions to the covering problem in selected dimensions. Values reported 
are taken from Conway and Sloane [la ], except for d = 7,8 and 9, which were obtained from 
Schiirmann and Vallentin [45j]. It is only in one and two dimensions that these solutions have been 
proved to be globally optimal [181 ] . 



Dimension, d 


Covering 


Covering Density, 9 


i 


A * 17 

A 1 = £t 


1 


2 


A * A 

A 2 = A 2 


1.2092 


n 
O 


^3 = ^3 


1 A CO CC 


4 


A * 

A 4 


1. moo 


e 



A * 


19/19 


6 


Li 


2.4648 


7 


L c 7 


2.9000 


8 


^8 


3.1422 


9 


A 5 


4.3401 


10 


A* 
^10 


5.2517 


12 


A* 


7.5101 


16 


A* 


15.3109 


17 


A* 


18.2878 


18 


A* 


21.8409 


24 


A 2 4 = A^ 


7.9035 



The best lower bound on the covering density is given by 

0>r d . 



(43) 



To define Td, let S be a regular simplex with edge length equal to two. Spheres of radius 
\f2dj (d + 1) centered at the vertices of S just cover S. The quantity is the ratio of the 
sums of the intersections of these spheres with S to the volume of S. Thus, in the large-c? 
limit, r — >■ d/e 3 ^ 2 . 



D. Quantizer Problem 

Consider a point process V in M d with configuration r 1; r 2 , . . . , r^. A <i-dimensional 
quantizer is device that takes as an input a point at position x in M. d generated from some 
probability density function p(x) and outputs the nearest point of the point process to x. 
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FIG. 4: (Color online) For lattices or periodic point process in which each point is equivalent, the 
Voronoi cells are congruent to one another, the Voronoi centroids coincide with the points of the 
point process, and the covering radius 1Z C is equal to the circumradius of the associated Voronoi 
cell. For an irregular point pattern, generally, the Voronoi cells are not congruent to one another, 
the Voronoi centroids do not coincide with the points of the point process, and and the covering 
radius 1Z C is not equal to the circumradius of the associated Voronoi cell. These two instances are 
illustrated in two dimensions for the triangular lattice and an irregular point pattern. 



• 


T • r t 


o 











FIG. 5: (Color online) Examples of two-dimensional quantizers. Any point x is quantized 
( "rounded-off" ) to the nearest point r-j. Left panel: Triangular lattice (best quantizer in M? fl8|). 
Right panel: Irregular point process. 



Equivalently, if the input x belongs to the Voronoi cell V(r-j), the output is r-j (see Fig. [5]) 
For simplicity, we assume that x is uniformly distributed over a large ball in M d containing 
the N points of the point process. One attempts to choose the configuration ri, T2, . . . , rjv 
of the point process to minimize the mean squared error, i.e., the expected value of |x — r^ 2 . 
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Specifically, the quantizer problem is to choose the Appoint configuration so as to minimize 
the scaled dimensionless error (sometimes called the distortion) 18] 



(44) 



where 



(R 2 



lim 

N^oo N 



1 N 

-Y 



8=1 



V(r<) 



(V(V)> 



1+: 



is a dimensionless error, 



(^(V)> 



1 * 



1=1 



(45) 



(46) 



is the expected volume of a Voronoi cell, and V(V(rj)) is the volume of the ith Voronoi 
cell. The scaling factor 1/d is included to compare the second moments appropriately across 
dimensions. We will denote the minimal scaled dimensionless error by Q m in- Note that since 
there is one point per Voronoi cell, the number density p is set equal to unity. 

If all of the Voronoi cells are congruent (as they are in the case of all lattices and some 
periodic point processes), we have the simpler expression 



(47) 



where the centroid of the Voronoi cell V is the origin of the coordinate system. The lattice 
quantizer 'problem is to find the lattice for which Q, given by (I47p . is minimum. Thus, Q can 
be interpreted as the scaled, dimensionless second moment of inertia of the Voronoi cell. 

The best known quantizers in any dimension d are usually lattices that are the duals 
of the densest known packings (see the discussion in Sec. IIII Al and Ref. 18]), except in 
dimensions 9 and 10, where the best solutions are still lattices, but are not duals of the 
densest lattice packings in those dimensions. Although the best known solutions of the 
quantizer and covering problems are the same in the first three space dimensions, they are 
generally different for d > 4 



451 ] . Zador [46j has derived upper and lower bounds on 
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Qmin- Conway and Sloane [18] have obtained conjectural lower bounds on Q m %n- We defer 
the discussion of these bounds to Sec. IVHt where we derive sharper upper bounds on Q m i n , 
among other results. Table IIVI provides the best known solutions to the quantizer problem 
in selected dimensions. 
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TABLE IV: Best known solutions to the quantizer problem in selected dimensions. Values reported 
are taken from Conway and Sloane 



and Eriksson 



la ], except for d = 9 and 10, which were obtained from Agrell 



471 ] . It is only in one and two dimensions that these solutions have been proved to 



be globally optimal 



18] 



Dimension, d 


Quantizer 


Scaled Error, Q 


i 


A * 17 

A 1 = £s 


0.08oodd 


2 


A * A 

A 2 = A 2 


0.080188 


Q 

o 


A* — n* 


U.U ( 0041:0 


4 


D\ = D A 


0.076603 


5 


D 5 * 


0.075625 


6 


E* 


0.074244 


7 


E* 


0.073116 


8 


Eg = Eg 


0.071682 


9 


T AE 


0.071626 


10 


Dto 


0.070814 


12 


K* 2 = K 12 


0.070100 


16 


AJ 6 = Aie 


0.068299 


24 


A24 = A24 


0.065771 



E. Comparison of the Four Problems 



Table |V] lists the best known solutions of the quantizer, covering and number- variance 
problems and sphere-packing problems in R d for selected d. It is seen that in the first two 
space dimensions, the best known solutions for each of these four problems are identical to 
one another. For d = 3, the densest sphere packing is the D3 or, equivalently, A3 lattice, 
which is the dual lattice associated with the best known solutions to the quantizer, covering 
and number-variance problems, which is the A% lattice. Thus, for the first three space 
dimensions, the best known solutions for each of the four problems are lattices, and either 
they are identical to one another or are duals of one another. However, such relationships 
may or may not exist for d > 4, depending on the peculiarities of the dimensions involved. 

There is a fundamental difference between the nature of the interactions for the sphere- 
packing problem and those for the other three problems. Any sphere packing (optimal 
or not) consisting of nonoverlapping spheres of diameter D is described by a short-ranged 
pair potential that is zero whenever the spheres do not overlap (when the pair separation 
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TABLE V: Comparison of the best known solutions to the quantizer, covering, number variance 
and sphere packing problems, as obtained from the previous four tables. Recall that the Es an d 
A24 lattices are self-dual lattices. 



Dimension, d 


Quantizer 


Covering 


Variance 


Packing 


1 


A\ = Z 


A\ = Z 


A\ = Z 


A\ = 7L 


2 


A* 2 = A 2 


A* 2 = A 2 


a; = a 2 


a; = a 2 


3 


Al = Dl 


Al EE D* 3 


Al EE Dl 


A 3 = D 3 


4 


D* 4 = D 4 


Al 


D* 4 = D 4 


Dl EE D A 


5 


D* 5 


A*, 


A? 




6 


E* 6 


L£ 


Et 


Eq 


7 


E* 7 


L, 7 


Af 


E 7 


8 


E 8 


Lt 


Es 


E 8 


9 


tAE 


Al 


Al 


A 9 


10 


E>to 


A* 


A* 


Pwc 


12 


K 12 


A* 


A max* 


A max 
A 12 


16 


Ke 


A* 


A* 


Aie 


24 


A 24 


A 2 4 


A 2 4 


A 24 



distance is greater than D) and is infinite whenever the pair separation distance is less than 
D. By contrast, the other three problems are described by "soft" bounded interactions. 
In particular, we have seen that the number variance is specified by a bounded repulsive 
pair potential with compact support [cf. (151]) ] . We will see in the subsequent sections that 
the covering and quantizer problems are described by many-particle bounded interactions 
but of the more general form (JT]), which involves single-body, two-body, three-body, and 
higher-body interactions. 

One simple reason why the optimal solutions of the sphere-packing and number-variance 
problems are related to one another either directly or via their dual solutions in the first three 
space dimensions is that they both involve short-ranged repulsive pair interactions only. The 
reader is referred to Refs. 23] and {5] for a comprehensive explanation. The reasons why the 
optimal solutions to these two problems are sometimes the optimal solutions for the covering 
and quantizer problems will become apparent in the subsequent sections. The explanation 
for why the optimal covering and quantizer solutions are generally different for d > 4 is 
discussed in Sec. |V] We note that the Leech lattice A24 for d = 24 is an exceptional case 
in that it provides the optimal solution to all four different problems. The remarkably high 
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degree of symmetry possessed by this self-dual lattice [22| accounts for this unique property. 
The only other dimensions where all four optimal solutions are the same are d — 1 and 
d = 2. 



IV. NEAREST-NEIGHBOR FUNCTIONS 



A. Definitions 



We recall the definition of the "void" nearest-neighbor probability density function Hy(R) 

m 

Hy(R)dR = Probability that a point of the point process lies at 

a distance between R and R + dR from a randomly (48) 
chosen point in ~R d . 

The "void" exclusion probability Ey(R) is the complementary cumulative distribution func- 
tion associated with Hy(R): 

POO 

E V (R) = / H v (x)dx ) (49) 
Jr 

and hence is a monotonically decreasing function of R Thus, Ey(R) has the following 
probabilistic interpretation: 

Ey(R) = Probability of finding a randomly placed spherical 

(50) 

cavity of radius R empty of any points. 

There is another interpretation of Ey that involves circumscribing spheres of radius R around 
each point in a realization of the point process. It immediately follows that Ey(R) is the 
expected fraction of space not covered by these circumscribing spheres. Differentiating f )49l 
with respect to R gives 

HrW = (51) 
Note that these void quantities are different from the "particle" nearest-neighbor functions 
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48 



49| in which the sphere of radius R is centered at an actual point of the point process 



(as opposed to an arbitrary point in the space). 
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It is useful to introduce the "conditional" nearest-neighbor function Gy(r) [30|, |3l|, which 
is defined in terms of Hy(r) and Ey(r) as follows: 



H V (R) = pSl (R)G v (R)E v (R), (52) 

where si(R) is the surface area of a ci-dimensional sphere of radius R [cf. (jSJ)]. Thus, we 
have the following interpretation of the conditional function: 

psi(R)Gv{R) dR = Given that a spherical cavity of radius R 

centered at an arbitrary point in the space is 
empty of any points of the point process, the 
probability of finding a point in the spheri- 
cal shell of volume Si(R) dR surrounding the 
arbitrary point. 

Therefore, it follows from (IST1) and (|52|) that the exclusion probability can be expressed in 
terms of Gy via the relation 

-R 



(53) 



Ey(R) = exp 



-psi(l) / x d ~ x G v {x) dx 



(54) 



It is clear that the void functions have the following behaviors at the origin for d > 2 [501 ]: 



Ev(0) = 1, H v (0) = 0, GV(0) = 1. (55) 

Moments of the nearest-neighbor function Hy(R) arise in rigorous bounds for transport 
properties of random media [2|. The nth moment of Hy(R) is defined as 

(R n )= R n H v {R)dR = n R n ~ l E v (R) dR. (56) 
Jo Jo 

B. Series Representations 

The void functions can be expressed as infinite series whose terms are integrals over the 



n-particle density functions |2j, |3j|. For example, the void exclusion probability functions 
for a translationally invariant point process are respectively given by 

00 k r ^ 

E v (R) = l + J2(-±) k lj / 9k(T U ...,T k )l[e(R-\x-T J \)dr j , (57) 
*=i fc - jRd j=i 



25 



where g n is the n-particle correlation function and Q(x) is the Heaviside step function defined 
by (|27p . The corresponding series for Hy(R) is obtained from the series above using flSTJ) . 
Note that the series (1571) can be rewritten in terms of intersection volumes of spheres: 

E V (R) = 1 + V(-l) fc ^r / <7 fc (r x , . . . , r fc ) vf(n, . . . , r fc ; R) dx x ■ ■ ■ dr k , (58) 



k=l 



where 



dx J] e(i?-|x- rj |) (59) 

is the intersection volume of n equal spheres of radius i? centered at positions ri, . . . , r n . 
Observe that v^fr; R)) = Vi(R)a(r; R), where V\{R) is the volume of a sphere of radius R 
[cf. ( 1141) ] and a(r;i?) is the scaled intersection volume [cf. (l2"2]l ] and fl24l)]. 

In the special case of a Poisson point distribution, g n = 1 for all n, and hence (|58|) 
immediately yields the well-known exact result for such a spatially uncorrelated point process 

E V (R) =exp(-pv 1 (R)). (60) 

The use of this relation with definition (|5T|) gives 

H V (R) = p Sl (R) expi-pv^R)). (61) 

For a single realization of N points within a large volume V in M d , we have 

E V {R) = l- pVl (R) + L^2 vf\r ij; R)-^J2 r ik , r jk - R) - ■ ■ ■ (62) 

i<j i<j<k 

This formula assumes that iV is sufficiently large so that boundary effects can be neglected. 
The second term in (16"2"|) pvi(R) = J2iLi v i(R)/V can be interpreted as a sum over one-body 
terms, which is independent of the point configuration. Clearly, the (n + l)th term in ( 16"2"|) 
can be interpreted as a sum over intrinsic n-body interactions,namely, t>™*. Thus, except 
for the trivial constant of unity (the first term), Ey(R) can be regarded to be a many-body 
potential of the general form ([!]), which heretofore was not observed. 



C. Rigorous Bounds on the Nearest-Neighbor Functions 

Upper and lower bounds on the so-called canonical n-point correlation function 
iJ n (x m ; x p_m ; r q ) (with n = p + q and m < p) for point processes in M. d have been found 
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0, |51|. Since the void exclusion probability and nearest-neighbor probability density func- 
tion are just special cases of H n , then we also have strict bounds on them for such models. 
Let X represent either Ey or Hy and represent the kth term of the series for these 
functions. Furthermore, let 

i 

W £ = (63) 

fc=0 

be the partial sum. Then it follows that for any of the exclusion probabilities or nearest- 
neighbor probability density functions, we have the bounds 

X < W\ for i even 

X > W\ for £ odd. (64) 

Application of the aforementioned inequalities yield the first three successive bounds on 
the nonnegative exclusion probability: 

E V (R) < 1 (65) 
E V {R) > l-pv 1 {R) (66) 

2 f2R 

E V (R) < 1 - pVl (R) + ^ Sl (l) J x d - 1 v$> t (x;R)g 2 (x)dx, (67) 

where v t { lt {x\ R) = vi(R)a(x; R) is the intersection volume of two d- dimensional spheres of 
radius R whose centers are separated by the distance x and a(x; R) is the scaled intersection 
volume given by fl22|) . The corresponding first two nontrivial bounds on the nonnegative 
pore-size density function H V (R) £1X6 cLS follows: 

H V {R) < p Sl {R) (68) 

2 r2R 

Hy{R) > p Sl (R) - / x d - 1 sf t (x;R)g 2 (x)dx, (69) 

* Jo 

where s 2 nt {x; R) = dv l 2 nt (x; R)/dR is the surface area of the intersection volume v l 2 nt {x\ R). 
Bounds on the conditional function Gy(r) follow by combining the bounds above on Ey{r) 
and Hy(r) and definition (1521) . For example, the following bounds have been found j^ : 

fMR) £ (70) 



and 



l- 7 ^ )Sl (l)Cx^sr^R)92(x) dx 
G v{R) > — . o — r2R - - — — - — — , (71) 



l-pv 1 (R) + £si(l) ^ R x d - 1 vi nt (x;R)g 2 (x) dx' 
which should only be applied for R such that Gy(R) remains positive. 



27 



D. Truncation of the Series Expansions for Nearest-Neighbor Functions for Pack- 
ings 



For congruent sphere packings of diameter D at packing density 0, the infinite series 
expansion for Ey(R) [cf. (j57j) ] will truncate after a finite number of terms for a bounded 
value of the radius R. A spherical region of radius R centered at an arbitrary point in the 
space exterior to the spheres can contain at most n max sphere centers. Therefore, series 
truncates after n max + 1 terms, i.e., 

E v (R) = l+J2(-V k Ti / gk(ri,... } T h )Y[Q{R-\x-T j \)dr j , (72) 

k=l ' J Rd j=l 

For a spherical region of radius R = D / a/3, n max = 2 and hence we have the exact 
expression that applies for < R < D / y/3 

2 r2R 

E V (R) = 1 - pv x {R) + ^-si(l) / x d - 1 vf t (x;R)g 2 (x)dx 

2 Jd 

= 1 - 2 d (£\ " + (J) " £ R x d -'a{x; R)g 2 (x)dx. (73) 

For any sphere packing for which fflTJl) applies such that r 2 /r 1 > 2/ \/3, we have upon use of 
(|73l) the exact result 

( 1 - 2 d <f) (%) d , 0<R<D/2 
E V (R) = I VKdJ , - - ' (74) 

y 1 _ 2 d <p [l - §a(D; R)] (§) d , D/2 < R > D/V3, 

where we have used the fact that r 1 = D. Importantly, this relation applies to the densest 
known lattice packings of spheres, at least for dimensions in the range 1 < d < 24. 



V. REFORMULATION OF THE COVERING AND QUANTIZER PROBLEMS 
A. Reformulations 

Now we can reformulate the covering and quantizers problems in terms of the void ex- 
clusion probability. In particular, the covering problem asks for the point process in M. d at 
unit density (p = 1) that minimizes the support of the radial function Ey(R). We define 
T^rmn ^ e smallest possible value of the covering radius 1Z C among all point processes for 
which Ey(R) = 0, which we call the minimal covering radius. This is indeed a special 
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ground state in which the "energy" is identically zero (i.e., Ey(lZ™ m ) = 0) |52j. Depending 
on the space dimension d, this special ground state will involve one-body interactions, the 
first two terms of ( 16 2 p . one- and two-body interactions, the first three terms of (IB"2"]) . one-, 
two- and three-body interactions, the first four terms of (16"2l . etc., and will truncate at 
some particular level, provided that Ey(R) for the point process has compact support. The 
minimal covering radius VJ^ m increases with the space dimension d and, generally speaking, 
the highest-order n-body interaction required to fully characterize the associated Ey(R) in- 
creases with d. Note that for a particular point process, twice the covering radius 2TZ C can be 
viewed as the "effective interaction range" between any pair of points, since the intersection 
volume t>2 nt (rjj; R), which appears in expression ( 16 2 p for Ey(R), is exactly zero for any pair 
separation > 2TZ C ; moreover, for such pair separations v 1 ^ can be written purely in terms 
of the lower-order intersection volume v^_ 1 . Therefore, because v 1 ^ > v 1 ^ for n > 3, the 
effective interaction range between any n points for n > 3 is still given by 2TZ C . 

The quantizer problem asks for the point process in M d at unit density that minimizes 
the scaled average squared error Q defined as 

g = (R 2 ) = / R 2 H v (R)dR = - / RE v (R)dR. (75) 
d d J d J 

We will call the minimal error Q m in- Thus, we seek the ground state of the many-body 
interactions that are involved upon substitution of fl6"2"l) into (1751) . Again, depending on 
the dimension, this many-body energy will truncate at some particular level, provided that 
Ey{R) has compact support. Again, as d increases, successively higher-order interactions in 
the expression (162]) must be incorporated to completely characterize Ey(R). 

B. Explicit Calculations for Some Common Lattices Using These Reformulations 

It is instructive to express explicitly the void exclusion probabilities for some common 
lattices and use these functions to evaluate explicitly their corresponding covering densities 
and scaled average squared errors in the first three space dimensions. 

In the simplest case of one dimension, the series expansion for Ey{R) [cf. (162]) ] for the 
integer lattice Z truncates after only one-body terms. At unit number density (p = 1), it is 
trivial to show that 

o, r > n c , 
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where the covering radius TZ C = 1/2, vi(R) = 2R, and the nearest-neighbor distance from 
a lattice point is unity. Using the definitions f )36|) and f!75|) in combination with f!76|) yield 
the covering density and scaled average squared error , respectively, for the optimal integer 
lattice: 

9 = 1, (77) 
Q = = 0.083333 .... (78) 



2 



Let us now determine the void exclusion probabilities for the Z (square) and A 2 = A 
(triangular) and lattices for R up to their respective covering radii for which Ey(R) = 0. 
For these lattices, the series expression ( |62l) for Ey(r) truncates after two-body terms. For 
the square and triangular lattices at p = 1, 



Ey(R) 



l-v^R), 0<R<n/2, 

1 - v^R) + 2v% A {r 1 ; R), n/2 <R<TZ C , (79) 

0, i? > ft c 



where ri is the nearest-neighbor distance from a lattice point and the covering radius 1Z C , 
equal to one half of the next-nearest-neighbor distance, which we will denote by r 2 . For 
the square and triangular lattices at p = 1, r x = 1 and 1Z C = y/2/2 = 0.7071..., and 
n = V2/3 1/4 = 1.0745 . . . and TZ C = V2/3 3/4 = 0.6204 . . . respectively. Figure |6] provides 
plots of Ey(R) for these two d = 2 lattices. Employing the definitions (13 6 j) and (l75p 
in combination with (I79|) provide the covering density and scaled average squared error, 
respectively, for the Z 2 lattice: 

= 1 = 1.57079..., (80) 

£ = -^ = 0.083333.... (81) 

Similarly, the corresponding equations for the A 2 lattice yields the covering density and 
scaled average squared error , respectively, for this optimal structure: 

2vr 



3V3 
5 



1.20919..., (82) 

_ = 0.08018.... (83) 
36v3 

It is also useful to express explicitly the void exclusion probabilities for the Z 3 (simple 
cubic) and A% (bcc) lattices for R up to their respective covering radii for which Ey(R) = 0. 
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FIG. 6: (Color online) The void exclusion probability Ey(R) for the Z 2 (square) and = A\ 
(triangular) lattice have support up to the covering radii 1Z C = v2/2 = 0.7071... and 1Z C = 
-v/2/3 3 / 4 = 0.6204 . . ., respectively, at unit number density (p = 1). 

It turns out that calculating Ey(R) in the case of the simple-cubic lattice is more complicated 
than that for the bcc lattice because the former involves up through four-body terms, i.e., 
t>4 nt . However, the symmetry of the geometry for Z 3 enables one to express v 1 ^ purely in 
terms of v 2 nt and v^. The calculation of Ey(R) does not involve four-body terms. For the 
simple cubic lattice at p = 1, 

' 1-vtiR), <R<n/2, 

1 - Vl (R) + 34 nt (r i; R), n/2 < R < r 2 /2, 

1 - Vl (R) + 3vf( ri ; R) + 34 nt (r 2 ; R) - 6^ nt (n, n, r 2 ; r 2 /2 < R < 1Z C , 

0, > n c , 

(84) 

where n = 1, r 2 = = 1.4142 . . ., K c = V3/2 = 0.8660 ... and vf{r, s, t; R) is explicitly 
given by (IBlj) in Appendix |B] for triangles of side lengths r, s and t. Here we have used 
the fact that for r 2 /2 < R < 1Z C , v^i^ri, r±, r 2 , r 2 ) = 6i>3 nt (ri, ri, r 2 ) — 3f 2 nt (r 2 ). Using the 
definitions (1361) and ( |75l) in combination with (1841) yield the covering density and scaled 
average squared error, respectively, for the Z 3 lattice: 

= ^5 = 2.72069..., (85) 



E V (R) 
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6 =12 



0.083333 .... 



(86) 



For the bcc lattice at p — 1, 



E V (R) 



l-vx{R), 0<R<n/2, 

1 - Vi(R) + 44 nt (r i; R), n/2 < R < r 2 /2, 

1 - + 44 nt (r i; i2) + 34 nt (r 2 ; r 2 /2 <R<R T , 

1 - + 44 nt (r i; R) + 3f 2 nt (r 2 ; i?) - 12u^*(ri,ri,r 2 ; R T < R < H B , 

0, i2 > ft c , 

(87) 

where n = v^/^/ 3 = 1.0911 . . ., r 2 = 4 1 / 6 = 1.2599 . . ., K c = v / 5/2 5 / 3 = 0.7043 . . ., and 
Rt = 3/2 13//6 = 0.6681 ... is the circumradius of a triangle of side lengths r±, r± and r 2 , the 
general expression of which is given by (1B2|) in Appendix [B] Employing the definitions fl36|) 
and (1751) in combination with (IHTjl provide the covering density and scaled mean squared 



error, respectively, for the A3 lattice: 



7T 



. 5 3/2 



24 
19 



1.46350... 



0.078543.... 



(88) 



192 ■ 2V3 

Figure [7| provides plots of Ey(R) for the simple cubic and bcc three-dimensional lattices 



(89) 



C. Remarks About Higher Dimensions 

We see that the best known solutions to the covering and quantizer problems are identical 
for the first three space dimensions. However, there is no reason to expect that the optimal 
solutions for these two problems to be the same in higher dimensions, except for d = 24 for 
reasons mentioned in Sec. IIII El Although both problems involve ground states associated 
with the "many-body interaction" function Ey(R), such that it possesses compact support 
for finite d, the precise shape of the function Ey(R) for a particular point configuration is 
crucial in determining its first moment or quantizer error. By contrast, the best covering 
seeks to find the point configuration that minimizes the support of Ey(R) without regard 
to its shape. 
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FIG. 7: (Color online) The void exclusion probability Ey(R) for the Z 3 (simple cubic) lattice and 
A% (bcc) lattice have support up to the covering radii K c = y/l/2 = 0.8660 . . . and K c = V5/2 5 / 3 = 
0.7043 . . ., respectively, at unit number density (p = 1). 

VI. RESULTS FOR THE COVERING PROBLEM 

Saturated sphere packings in M. d should provide relatively thin coverings. Surrounding 
every sphere of diameter D in any saturated packing of congruent spheres in M. d at packing 
density <j) s by spheres of radius D provides a covering of M. d , and thus the associated covering 
density 9 S is given by 

e s = PsVl (D) = 2%, (90) 

where p s and <p s = p s v\(D/2) are the number density and packing density, respectively, of 
the saturated packing. What is the thinnest possible covering associated with a saturated 
packing? It immediately follows that the thinnest coverings among saturated congruent 
sphere packings in M. d are given by the saturated packings that have the minimal packing 
density cj>* in that space dimension and have covering density 

0: = 2U:. (91) 

We can also bound the packing density <p s of any saturated sphere packing in ~R d from 
above using upper bounds on the covering density. 

Lemma 1: The density <ft s of any saturated sphere packing in M. d is bounded from above 
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lattice 


a; 


lattice 



according to 



5d rfln(d) dln(ln(d)) 
S 2 d + 2 d + 2 d ' 1 j 



The proof is trivial in light of the upper bound ( 142|) on the covering density and relation 
PJ). 

The standard random sequential addition (RSA) sphere packing is a time-dependent 



process produced by random 
into a large region V C M. d 



y, irreversibly, and sequentially placing nonoverlapping spheres 



531 ] . Initially, this large region is empty of sphere centers and 



subsequently spheres are added provided each attempted placement of a sphere does not 
overlap an existing sphere in the packing. If an attempted placement results in an overlap, 
further attempts are made until the sphere can be added to the packing at that time without 
violating the impenetrability constraint. For identical <i-dimensional RSA spheres, the filling 
process terminates at the saturation limit at infinitely long times in the infinite- volume limit. 
Thus, in this limit, the RSA packing is a saturated packing. The saturation density <p s can 
only be determined exactly in one dimension, where it is known to be <p s = 0.747597. . . 

. In higher dimensions, the saturation density can only be determined from computer 
simulations. Earlier work focused on two and three dimensions, where is was found that 
<p s ~ 0.547 55j| and <p s ~ 0.38 [56], respectively. More recent numerical work has reported 



RSA saturation densities for the first six space dimensions [38] . It is important to emphasize 
that upper bound (192]) provides a very poor estimate of the RSA saturation density; for 
example, for d = 3, it gives an upper bound greater than unity (2.32224, which is about 6.1 
times larger than the value obtained from simulations) and for d — 6, it provides the upper 
bound 0.691402, which is about 7.3 times larger than the value obtained from simulations 



38| . Thus, bound (192]) could only be realized by RSA saturated packings, if at all, in high 



dimensions. 

Therefore, it is useful here to compute the corresponding covering densities for RSA 
packings. The numerical data for RSA saturation densities for 2 < d < 6 reported in 
Ref. [38] was well approximated (with a correlation coefficient of 0.999) by the following 
form: 

#. = £ + g. <*» 
where C\ = 0.202048 and c<i = 0.973872. This can be used to estimate the RSA saturation 
densities for d > 7. However, even though the upper bound (192]) on the packing density 
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of a saturated sphere packing grossly overestimates the RSA value for the first six space 
dimensions, we include in the fit function for <p s a dhx(d) correction, which will lead to a 
more conservative estimate of the covering density, as explained shortly 57| . Including such 
a term, we find the following fitted function for the saturated RSA packing density: 

a x a 2 d a 3 d\n(d) 
h = & + 25" + & (94) 

with a correlation coefficient of 0.999, where a x = 0.350648, a 2 = 0.87660 and a 3 = 0.041428, 
does as well as ( )93l) for 2 < d < 6. Since (|94l) predicts slightly higher densities than (|93l) for 
7 < d < 24 58] , we use it to obtain the corresponding estimate of the RSA covering density, 
namely, 

9 S = ai + a 2 d + asdln(d), (95) 

which is a slightly more conservative estimate, since ( 193]) would yield a slightly thinner 
covering density. 

In Table IVIj we provide estimates of the RSA covering densities for selected dimensions up 
through d = 24. The covering density for d — 1 is determined from Rehyi's exact saturation 
packing density value [54j and multiplying it by 2. The values for 2 < d < 6 is are obtained 
from the reported saturated density values in Ref. [38( and multiplying each density value 
by 2 d . The values reported for d > 7 are estimates obtained from the fitting formula 
( |95i) . Comparing Table [VTl to Table [TTT1 for the best known coverings, we see that saturated 
RSA packings not only provide relatively thin coverings but putatively represent the first 
non-lattices that yield thinner coverings than the best known lattice coverings beginning in 
about dimension 17. This suggests that saturated RSA packings may be thinner than the 
previously best known coverings for 17 < d < 23 and probably for some dimensions greater 
than 24. 



VII. RESULTS FOR THE QUANTIZER PROBLEM 

Using the successive lower and upper bounds on the void exclusion probability function 
Ey(R) given in the previous section, we can, in principle, derive corresponding bounds on 
the minimal error Q m in- Moreover, one can obtain a variety of upper bounds on Q m %n using 
our approach by utilizing the exact form of the void exclusion probability, when it is known, 
for some point process at unit density. Since a general point process must have an error Q 
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TABLE VI: Covering density 8 S for RSA packings at the saturation state in selected dimensions. 
The values for the first 6 space dimensions are obtained from the reported saturated density values 
given in Refs. 



38, 



541 ] and using (|90p . The values reported for d > 7 are estimates obtained from the 



fitting formula (|95p . Included in the table are the corresponding the saturation packing densities. 



Dimension, d 


Covering Density, 6 S 


Packing Density, (f) s 


-i 

i 


1.4952 


0.74759 


2 


2.1880 


0.54700 


3 


o.Uozz 


f\ 00070 

O.o8z78 


A 

4 


4.0726 


0.25454 


c 



0. 10ZO 


U.101UZ 


6 


6.0121 


0.09394 


7 


7.0512 


0.05508 


8 


8.0526 


0.03145 


9 


10.0706 


0.01769 


10 


11.0860 


0.009834 


12 


12.1052 


0.002955 


16 


16.2141 


2.4740 x 10" 4 


17 


17.2482 


1.3159 x 10" 4 


18 


18.2848 


6.9751 x 10~ 5 


24 


24.5489 


1.4632 x lO" 6 



that is generally larger than the minimal Q m ini it trivially follows that 

Qmin — Q • 



(96) 



A. Revisiting Zador's Bounds 



To illustrate how we can obtain bounds on Q m i n using our approach, we begin by red- 
eriving the following bounds due to Zador |46| : 



— i— r(i + d/2) 2 > d < g mm < J-r(i + d/2) 2 / d v(i + 2/d). 

[a + 2)7r an 



(97) 



Consider the lower bound first. Combination of relation (J75J) and lower bound (|66p yields 
at unit density 



2 r R ° 1 

Qmin > 1 R[l~ V X {R)]dR = — — -r(l + d/2fl\ 

d Jo [d + 2)7T 



(98) 
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which is seen to be equal to Zador's lower bound. Here Rq = T(l + d/2) 1 / d /y/n is the 
zero of 1 — V\(R). It is clear that a sphere of radius R$ has the smallest second moment of 
inertia of any solid d- dimensional solid, and hence establishes the lower bound. The simplest 
example of a point process for which Ey(R) is known is the Poisson point process [cf. (160))] . 
Substitution of ( 160]) into ( 196]) at unit density yields 

2 r°° 1 

Qmva < ^ / Rexpi-v^dR = — r(l + d/2) 2 ' d T(l + 2/d), (99) 
a J dir 

which is seen to be equal to Zador's upper bound. These derivations of the inequalities 
stated in ( 1971) appear to be much simpler than the ones presented by Zador. 

In the lenge-d limit, Zador's upper and lower bounds become identical, and hence one 
obtains the exact asymptotic result 

G min -> — = 0.058550 ... as d oo. (100) 
2ne 

The convergence of Zador's bounds to the exact asymptotic limit is to be contrasted with the 
sphere-packing problem in which the best upper and lower bounds on the maximal density 
become exponentially far apart in the high-dimensional limit. 



B. Improved Upper Bounds 

Improved upper bounds on Q m %n can be obtained by considering those point processes 
corresponding to a sphere packing for which the minimal pair separation is D and lower 
bounds on the conditional function Gv{R) for R > D/2. In what follows, we present two 
different upper bounds on Q m in based on this idea that improve upon Zador's upper bound. 

For any packing of identical spheres with diameter D, the following exact relations on 
the nearest-neighbor quantities apply for R < D/2 

M*> = !-*#(§)'. ffvW^g)", G ^) = T -^ |f 0<R<D,2. 

(101) 

Observe that Gv{R) is a monotonically increasing function of R in the interval [0, D/2]. 
From these equalities, it immediately follows that 

Evip/2) = !-<!>, H v (D/2) = ^, G v (D/2) = -^—. (102) 
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Consider now the class of sphere packings for which the conditional nearest-neighbor 
function is bounded from below according to 

G V (R) > — — for all R > D/2. (103) 
1-0 

This class of packings includes equilibrium (Gibbs) ensembles of hard spheres along the 



disordered fluid branch of the phase diagram 0, [|o], 48j , nonequilibrium disordered sphere 
packings, such as the "ghost" random sequential addition (RSA) process jig , and a large 
class of lattice packings of spheres, as will be described below. In the equilibrium cases, it is 
known that Gv{R) is a monotonically increasing function of R for d > 2 and thus using this 
property together with the equality Gy(D/2) = 1/(1 — 0) [cf. fl 1 2 [) ] means that the lower 
bound (11031) is obeyed. For one-dimensional equilibrium "rods," the bound (11031) is sharp 
(exact) Gy(R) for all realizable G [0, 1]. A bound of the type (11031) was used to bound the 
related "particle" mean nearest-neighbor distance from above for different classes of sphere 



packings for all d 48]. 



Using definition fl54l) and inequality (I103p . the void exclusion probability function obeys 
the following upper bound for R > D/2: 



E V (R) < (l-0)exp 



1 



R\ d 1 



D 2 d 



R>D/2. (104) 



Since any upper bound on the nonnegative function E V (R) leads to an upper bound on its 
first moment, we then have upon use of (196]) . (llOip and (I104p . the upper bound 



, 4j0T£+d/2]]^ 
dir 



(d + 2(1-0)) , (1-0) Sl- 



id 



exp 



4(2 + d) 2d V <P ) V 1 -^/ \d'\- 

(105) 

where T (s, x) = t s ~ 1 e~ t dt is the incomplete gamma function. Observe that the prefactor 
multiplying the bracketed expression is D 2 /d, where, in light of f JT4^) and ffl6|) . D = 2[0r(l + 
d/2)] 1 / d / y/n, assuming unit number density. Note also that the upper bound fl 105[) depends 
on a single parameter, namely, the packing density <p. Thus, there is an optimal packing 
density (p opt G [0, ma x] that yields the best (smallest) upper bound for any particular d, 
where 4> ma x is the maximal packing density. Since the right side of the inequality is a 
monotonically decreasing function of <fi for any d, then the optimal density opt is, in principle, 
given by (p max . It is noteworthy that the upper bound fl 105[) for the optimal choice = 4> ma x 
may still be valid for a packing even if the bound (I104p . upon which it is based, is violated 
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for R of the order of D because the exponential tail can more than compensate for such a 
violation such that the error [first moment of Ey(R)] is overestimated. Observe also that 
because 

/9 rh \ 1 T f\ - rh\ 2/d 

+ 0(1) (rf^oo) (106) 



1 

d 1 l-<j)J ~ 4 

the upper bound (11051) tends to the exact asymptotic result (11001) of (20e) _1 . 

Before discussing the optimal bounds, it is useful to begin with an application of the 
upper bound (11051) for the sub-optimal case of a disordered sphere packing, namely, the 
aforementioned ghost RSA packing process 15] , which we now show generally improves on 
Zador's upper bound. This represents the only exactly solvable disordered sphere-packing 
model for all realizable densities and in any dimension, as we now briefly describe. The 
ghost RSA packing process involves a (time-dependent) sequential addition of spheres in 
space subject to the nonoverlap condition. Not only is an attempted addition of a sphere 
rejected if it overlaps an existing sphere of the packing, it is also rejected if it overlaps any 
previously rejected sphere (called a "ghost" sphere). Unlike the standard RSA packing, 
the ghost RSA packing does not become a saturated packing in the infinite-time limit. All 
of the n-particle correlation functions for this nonequilibrium model have been obtained 
analytically for any d, time t, and for all realizable densities. For example, one can show 
that the maximal density (achieved at infinite time) is given by 

<t> = 5? ( 107 ) 

and the associated pair correlation function is 

, , N 26(r — D) , , 

"* (r ) = J^wy (108) 

where Q(x) is the unit step function, equal to zero for x < and unity for x > 1. It is 
straightforward to verify that the upper bound on the exclusion probability Ey(R) for this 
infinite-time case obtained by using fl 108[) in the inequality (157)1 is always below the upper 
bound ( 1104p . Therefore, the upper bound (11051) is valid at the maximal density, i.e., at 
cf) = l/2 d , we have 

[T(l + d/2)} 2 / d 



< 



dix 



(d + 2(l - l/2 d )) 2(1 - i/2 rf )( 2 + d )/ d , . x r 



4(2 + d) d * V2 d -ly \d'2 d -lj 

(109) 
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For d = 1, 2 and 3, this upper bound yields 0.166666..., 0.124339... and 0.106797..., 
respectively, which is to be compared to Zador's upper bound, which gives 0.5, 0.159154 . . . 
and 0.115802..., respectively. We note that in the large-d limit, the upper bound (I109p 
yields the exact asymptotic result f llOOp . which implies that the upper bound ( I104p on 
E V (R) becomes exact for ghost RSA packings, tending to the unit step function in this 
asymptotic limit, i.e., 

E v (R)^Q(r-D) (d-»oo). (110) 

This asymptotic result implies the following corresponding one for the void nearest-neighbor 
probability density function: 

H v (R)^5(r-D) (d->oo). (Ill) 
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FIG. 8: (Color online) Comparison of the exact result for Ey scaled by 1 — (p m ax for the optimal 
lattice packings for d = 2 (A* 2 ) and d = 4 (Z> 4 ), as obtained from flUD, for 1/2 < R/D < l/\/3 
(solid curves) to the corresponding estimates obtained from (|104p for these cases (dashed curves). 
It is only for the case d = 2 that estimate (|104p is not a rigorous pointwise upper bound on the 
exact void exclusion probability for 1 < d < 24 and, likely, for d > 24. The exponential tail 
associated with ()104|) more than compensates for the narrow pointwise violation for the special 
case d = 2, resulting in a strict upper bound on the first moment of Ey(R), i.e., (|105p remains a 
strict upper bound for d=2. 
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We now return to finding the optimal (smallest) upper bound (11051) for each dimension. 
For d — 1, the optimal density is (p op t = 4> ma x = 1, which produces the sharp bound 

Qmin < jo = 0.083333 ... . (112) 

This bound is exact in this case because the inequality f 1 1 3 j) is exact for all realizable 
densities for equilibrium hard "rods," including at <fi = 1, which corresponds to the optimal 
integer lattice packing. This is to be contrasted with Zador's upper bound, which yields 
1/2 for d = 1 and is far from the exact result. The improved upper bound f 1 1 5 p in the 
higher dimensions reported in Table IVHI is obtained by evaluating it at the densities of 
the densest known lattice packings in these respective dimensions 18|. We note that it is 
only for optimal triangular lattice packing in M. 2 that the upper bound (11041) on Ey(R) is 
violated pointwise for a small range of R around R/D = 1/2 [inequality fll04j) is obeyed for 
R/D > 0.539 and in the vicinity of R/D = 0.5], but the exponential tail associated with 
(11041) more than compensates for this narrow pointwise violation, resulting in a strict upper 
bound on the first moment of E V (R), i.e., (I105p remains a strict upper bound for d = 2. 
Using relation f[T4"|) and lattice coordination properties, it is easily verified that the (I104p is a 
strict upper bound on the void exclusion probability for the densest known lattice packings 
for all R > D/2 and 3 < d < 24 as well as d = 1, and hence inequality fl 105[) provides a 
strict upper bound on the scaled error for all of these lattices. For illustration purposes, 
we compare in Fig. [8] the exact result for Ey obtained from (1741) to the estimate (I104p for 
the cases d = 2 and d = 4 for 1/2 < R/D < 1/^3. The upper bound (11051) is generally 
appreciably tighter than Zador's upper bound for low to moderately high dimensions. 

C. Results for Saturated Packings 

For any saturated packing of identical spheres of diameter D, Ey(R) by definition is 
exactly zero for R beyond the diameter, i.e., 

E V (R) = R > D. (113) 

In the particular case of saturated RSA packings, the void exclusion probability can com- 



puted using the same techniques described in Ref. [38| for the first six space dimensions. 
These results are summarized in Fig. [91 The corresponding quantizer errors for these di- 
mensions are listed in Table IVIII1 We see that the discrepancies between the saturated 
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TABLE VII: Comparison of the best known quantizers in selected dimensions to the conjectured 
lower bound due to Conway and Sloane and the improved upper bound (|105p . 



d 


Quantizer 


Scaled Error, Q 


Conjectured 
Lower bound 


Improved 
Upper Bound 


1 


A\ = Z 


0.083333 


0.083333 


0.083333 


2 


A* = A 2 


0.080188 


0.080188 


0.080267 


3 


A* EE D* 


0.078543 


0.077875 


0.079724 


4 


D\ EE D 4 


0.076603 


0.07609 


0.078823 


5 


A|* 


0.075625 


0.07465 


0.078731 


6 


E* 


0.074244 


0.07347 


0.077779 


7 


Af 


0.073116 


0.07248 


0.076858 


8 


Eg = E 8 


0.071682 


0.07163 


0.075654 


9 


t AE 


0.071626 


0.070902 


0.075552 


10 


Dto 


0.070814 


0.070405 


0.074856 


12 


K* 12 = K 12 


0.070100 


0.06918 


0.073185 


16 


A*6 = A 16 


0.068299 


0.06759 


0.070399 


24 


^24 = ^24 


0.065771 


0.06561 


0.067209 




R/D 

FIG. 9: (Color online) The void exclusion probability Ey(R) for saturated RSA packings of 
congruent spheres of diameter D for the first six space dimensions. 

RSA quantizer error improves as d increases as compared to the best known quantizer error 
reported in Table HVl 

Lemma 2: Saturated sphere packings in M. d possess void nearest-neighbor functions that 
tend to the following high-dimensional asymptotic behaviors: 

Ev(R) -> 0(r - D), H v (R)^8(r-D) (d-»oo). (114) 
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TABLE VIII: The quantizer errors for saturated RSA packings in the first six space dimensions. 
We denote by Q s the quantizer error for such as saturated packing. 



Dimension, d 


Quantizer Error, Q s 


1 


0.11558 


2 


0.09900 


3 


0.09232 


4 


0.08410 


5 


0.07960 


6 


0.07799 



For any saturated packing at packing density <f> 8 , it is clear that Ey(R) is bounded from 
above for R > D/2 as follows: 



E V (R) < 1 - (f) s , D/2<R<D. 



(115) 



Let Q s denote the scaled dimensionless quantizer error for a saturated packing. Combination 
of the expression for E V (R) in (1101)) and (jl 15)) yields the following upper bound on Q s : 

4^r(i + d/2)] 2 / d 



dn 



(d+2(l- <(>.)) 3(1 



(116) 



4(2 + d) 8 

The fact that this upper bound becomes exact in the high- dimensional limit (that is, it tends 
to (27re) _1 ), implies that Ey(R) and Hy(R) for a saturated packing tends to the unit step 
function and radial delta function, as specified by ( 1114)) . Not surprisingly, the bound (1116)) 
is not that tight in relatively low dimensions. 



VIII. CONCLUDING REMARKS AND DISCUSSION 

We have reformulated the covering and quantizer problems as the determination of the 
ground states of interacting point particles in M. d that generally involve single-body, two- 
body, three-body, and higher-body interactions; see Sec. |V] The n-body interaction is di- 
rectly related to a purely geometrical problem, namely, the intersection volume of n spheres 
of radius R centered at n arbitrary points of the system in ~R d . This was done by linking 
the covering and quantizer problems to certain optimization problems involving the "void" 
nearest-neighbor functions. This reformulation allows one to employ theoretical and numer- 
ical optimization techniques to solve these energy minimization problems. 
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A key finding is that disordered saturated sphere packings provide relatively thin coverings 
and may yield thinner coverings than the best known lattice coverings in sufficiently large 
dimensions. In the case of the quantizer problem, we derived improved upper bounds on the 
quantizer error that utilize sphere-packing solutions. These improved bounds are generally 
substantially sharper than an existing upper bound due to Zador in low to moderately large 
dimensions. Moreover, we showed that disordered saturated sphere packings yield relatively 
good quantizers. Our reformulation helps to explain why the known solutions of quantizer 
and covering problems are identical in the first three space dimensions and why they can 
be different for d > 4. In the first three space dimensions, the best known solutions of the 
sphere packing and number variance problems are directly related to those of the covering 
and quantizer problems, but such relationships may or may not exist for d > 4, depending on 
the peculiarities of the dimensions involved. It is clear that as d becomes large, the quantizer 
problem becomes the easiest to solve among the four ground-state problems considered in 
this paper, since, unlike the other three problems, the asymptotic quantizer error tends to 
the same limit independent of the configuration of the point process. 

The detection of gravitational waves from various astrophysical sources has and will be 



searched for in the output of interferometric networks [59| by correlating the noisy output 
of each interferometer with a set of theoretical waveform templates [26| . Depending upon 
the source, the parameter space is generally multidimenisonal and can be as large as d — 17 
or larger for inspiraling binary black holes. The templates must cover the space and the 
challenge is to place them in some optimal fashion such that the fewest templates are used 
to reduce computational cost without reducing the detectability of the signals. 

Optimal template placement for gravitational wave data analysis has proved to be highly 
nontrivial. One solution proposed for the optimal placement in flat (Euclidean) space is to 



simply use the optimal solution to the covering problem |26[. However, this requires every 
point in the parameter space to be covered by at least one template, which rapidly becomes 
inefficient in higher dimensions when optimal lattice covering solutions are employed. An- 
other approach to template placement consists in relaxing the strict requirement of complete 
coverage for a given mismatch, and instead require coverage only with a certain confidence 



26l |. Such stochastic approaches have involved randomly placing templates in the parameter 



space, accompanied by a "pruning" step in which "redundant" templates, which are deemed 
to lie too close to each other, are removed 60|. The pruning step may be a complication 
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that can be avoided, as discussed below. Another approach is to place spherical templates 



down according to a Poisson point process 



611 ]. which has been claimed to provide good 



solutions for d > 10. The problem with the latter approach is that there will be numerous 
multiple overlaps of templates, which only increases as the number density of templates 
(intensity of the Poisson point process) is increased in order to cover as much of the space 
as is computationally feasible. 

The results of the present study suggest alternative solutions to the optimal template 
construction problem. First, we remark that if the covering of space by the templates is 
relaxed, then it is possible that the optimal lattice quantizers could serve as good solutions 
in relatively low dimensions (d < 10) because the mean square error is minimized. Second, 
in such relatively low dimensions, we have shown that saturated RSA sphere packings pro- 
vide both relatively good coverings and quantizers, and hence may be useful template-based 
constructions for the search of gravitational waves. Indeed, we have shown that saturated 
RSA sphere packings are expected to become better solutions as d becomes large. However, 
for d > 10, it will be computationally costly to create truly saturated packings, which by 
definition provide coverings of space (see Sec. IVII) . However, the existing stochastic ap- 
proaches do not require complete coverage of space and hence an unsaturated RSA sphere 
packing that gets relatively close to the saturation state might still be more computation- 
ally efficient than either the random placement /pruning technique 60j (because pruning is 
unnecessary) or the Poisson placement procedure 6l| (because far fewer spheres need to be 
added). Moreover, when the template parameter space is curved, which occurs in practice 
6l| , the RSA sphere packing would be computational faster to adapt than lattice solutions. 

In future work, we will explore whether our reformulations of the covering and quantizer 
problems as ground-state problems of many-body interactions of the form (CQ) can facilitate 
the search for better solutions to these optimization tasks. Clearly, a computational challenge 
in high dimensions will be the determination of the intersection volume v 1 ^ of n spheres of 
radius R at n different locations in M. d for sufficiently large n. However, it is possible that 
that the series representations ( 162|) for Ey(R) and bounds on this quantity [cf. Sec. IIV C] 
can be used to devise useful approximations of the monotonic function Ey(R), which should 
be zero for R > 27Z C , where 1Z C is the bounded covering radius. Such approximation could 
be employed to evolve an initial guess of the configuration of points within a fundamental 
cell to useful but sub-optimal solutions, which upon further refinement could suggest novel 
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solutions. In short, the implications of our reformulations to discover better solutions to the 
covering and quantizer problems in selected dimensions have yet to be fully investigated and 
deserves future attention. 
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Appendix A: Computing the Number Variance via the Epstein Zeta Function for 

d = 12, 16 and 24 

Here we summarize the steps in computing the asymptotic number-variance coefficient 
fl34|) for the lattices Ki%, Axe and A24, which correspond to the densest lattice packings in 
dimensions 12, 16 and 24, respectively. Importantly, the sum in (|34|) converges slowly. We 
noted in Sec. IIII Bl that the dual oi the lattice that minimizes the Epstein zeta function Z\(s) 
[defined by (1351) ] at s — (d + l)/2 among all lattices will minimize the asymptotic number- 
variance coefficient ( }34"1) among lattices. We will exploit number-theoretic representations 
of the Epstein zeta function that enable its efficient numerical evaluation and thus efficient 
computation of the asymptotic number-variance coefficient (134]) using the aforementioned 
duality relation. 

First, let us note that Epstein zeta function Z^{s) (|35|) can be rewritten as follows: 



where Zi is the coordination number at a radial distance Pi from some lattice point in the 
lattice A. [Note that the Epstein zeta function defined in this way applies to a general 
periodic point process provided that Zi is interpreted in the generalized sense discussed 




(Al) 
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in (jni)-] The quantities Zi and pi for many well-known lattices in M. d can be obtained 
analytically using the theta series for a lattice A, which is defined by 

oo 

e A (p) = i + 5> <g »*, (A2) 



1=1 



and is directly related to the quadratic form associated with the lattice 



. This series 



18| 



expression can usually be generated from the simpler functions 8 2 , 3 , and 84, which are 
defined by Q: 

=2Er ? (m+1/2)2 (A3) 
0M =1 + ^1^ (A4) 

eM =i + 2Et=i(-i) m2 - (as) 

Specifically, for the K 12 , A 16 and A 2 4 lattices, the associated theta series are given by 

Q Kl2 = 0o(2g) 6 + 450 o (2g) 2 0i(2g) 4 + 180 1 (2g) 6 

= l + 756g 4 + 4032g 6 + 20412g 8 + --- (A6) 

9 Al6 = ^[9 2 (2q) 1G + 9 3 (2q) ie + e,(2q) 1(i + 309 2 (2q)%(2q) s ] 

= 1 + 4320g 4 + 61440g 6 + • • • (A7) 
1 45 

e A24 = -Mg) 8 + e 3 (q) 8 + ^(<?) 8 ] 3 - -Mq)8M0Mf 

= 1 + 196560g 4 + 16773120g 6 + • • • (A8) 

where 

0o = 9 2 (2q)9 2 (6q) + 9 3 (2q)9 3 (6q), (A9) 

0! = 9 2 (2q)9 3 (6q) + 9 3 (2q)9 2 (6q). (A10) 

Direct evaluation of flAll) has the same convergence problems that the direct evaluation 
of the asymptotic number- variance coefficient fl34|) . However, we can exploit alternative 
number-theoretic representations of (1A1[) to facilitate its evaluation. In particular, there is 
an expression for the Epstein zeta function that can be derived using Poisson summation 

F A (s) = n- s T(s)Z A (s) 

= ^7--+ E G ( S ^IP| 2 )+ E G(^- S ,vr|p| 2 ), (All) 

S 7T S ^ 

2 P e A p g A* 

Pt^O p^O 
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where 

G(8,x) = x-'T(8,x) (A12) 

and T(s,x) = J x °° e~ t t s ~ 1 dt is the complementary incomplete Gamma function. It is impor- 
tant to note that the volumes of the fundamental cells of the lattice and its dual associated 
with the first and second sums in ( 1A11I) . respectively, are both taken to be unity here. Us- 
ing the appropriate theta series given above for the lattices corresponding to the densest 
lattice packings for d = 12, 16 and 24 and expression (lAllj) for s = (d + l)/2, one finds the 
corresponding Epstein zeta functions to be 

Z Kl2 (s = 13/2) = 12.527470092112 . . . (A13) 
Z Al6 (s = 17/2) = 2.606378060701 . . . (A14) 
Z A24 (s = 25/2) = 0.026464258871 . . . (A15) 

Now since all of these lattices are self-dual (i.e., Kyi = K{ 2) Ai 6 = A* 6 , A 2 4 = A24), we can 
directly determine the corresponding asymptotic number-variance from (|34|) by replacing 
the sum therein with the appropriate evaluation of the Epstein zeta function specified by re- 
lations (1A13j) - (lA15j) . The asymptotic number-variance values for d = 12, 16 and 24 reported 



in Table [TT] were obtained in this fashion. 



Appendix B: Intersection Volume of Three Spheres in Three Dimensions 



For d — 3, the intersection volume t>3 nt (r, s, t; R) of three identical spheres of radius R 



whose centers are separated by the distances r, s, and t for R > R? is given by 



62) 



4>,M;fl) = |-t+|tan-(- 2 + s2 + (2 _ 8JJ2 

- s (i?- s 7l2)tan-> ( r2 _ 2 y +j2 
-((tf-iV^taiW 2QrS 



Q ■ rst 



>H + s — 



where < tan 1 x < it, 
Rx = 



rst 



[(r + s + t)(-r + s + t)(r - s + t)(r + s - t)] 1 / 2 



(Bl) 



(B2) 
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is the circumradius of the triangle with side length lengths r, s, and t and Q = 
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